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We have investigated the electronic and magnetic response ol a single and a pair of interacting Fe 
atoms placed in patterned dehydrogenated channels in graphane within the framework of density 
functional theory. We have considered two channels - "armchair" and "zigzag" channels. Fully 
relaxed calculations have been carried out for three different widths of the channels. Our calculations 
reveal that the response to the magnetic impurities is very different for these two channels. We have 
also shown that one can stabilize magnetic impurities (Fe in the present case) along the channels 
of bare carbon atoms, giving rise to a magnetic insulator or a spin gapless semiconductor. Our 
calculations with spin-orbit coupling shows a large in-plane magnetic anisotropy energy for the case 
of an armchair channel. The magnetic exchange coupling between two Fe atoms placed in the 
semiconducting channel with an armchair edge is very weakly ferromagnetic whereas a fairly strong 
ferromagnetic coupling is observed for reasonable separations between Fe atoms in the zigzag edged 
metallic channel with the coupling mediated by the bare carbon atoms. The possibility of realizing 
an ultra-thin device with novel magnetic properties is discussed. 

PACS numbers: 71.15.Mb, 73.22.Pr, 73.20.At, 81.05.U- 



I. INTRODUCTION 

In the last few years graphene has generated an intense 
amount of activities due to its numerous unique proper- 
ties. [U-il Apart from a very strong multidisciplinary 
interest in the properties of pure graphene, the scien- 
tific community is engaged heavily in the modifications 
of graphene to realize new and novel effects. There is a 
tremendous interest in tuning the electronic structure of 
graphene to open up an energy gap suitable for the ap- 
plications in electronic devices. One of the ways is by the 
chemical modification of graphene [1, Q and the other is 
to realize confined structures in the form of nanoribbons, 
nanodots etc. [B-Q In this context, a considerable re- 
search effort is continuing in attaching hydrogen [Tol - [l5j 
or fluorine atoms [T6l42]| to graphene in order to have a 
stable composite system with a desired energy gap. 

Graphane is a hydrogenated graphene lattice with one 
hydrogen atom attached to each carbon atom giving rise 
to an insulating system with sp'^ bonds, first predicted 
by ab initio theory Q and then confirmed by experi- 
ments. The most attractive aspect of hydrogenation 
in graphene is the opening of a band gap. In a previous 
paper, we studied various concentrations of hydrogena- 
tion on graphene. We showed that a semimetal to metal 
to insulator transition occurs as a function of hydrogen 
concentration in graphene and hydrogenation favors clus- 
tered coirfigurations leading to the formation of compact 
islands. [22| 

One of the results in our earlier work as well as work 
reported by others [22| - [23 l showed that the patterning 
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of graphane via hydrogen desorption, the partial hydro- 
genation with the simultaneous presence of bare and hy- 
drogenated carbon atoms, is a novel way of modifying 
the properties of pure graphene. It turns out that such a 
pattern can lead to conducting channels, quantum dots, 
opening of band gap and finally magnetically coupled 
graphene/ graphane interfaces. These interesting proper- 
ties of graphene/graphane interface structures have at- 
tracted a large number of researchers. [23l - l27j Recent 
studies show that the energy band gaps of both arm- 
chair and zigzag graphane nanoribbons decrease as the 
nanoribbons become wider. [2^ A variety of other 
interesting properties such as thermal stability, interface 
engineering, thermal conductivity, etc. have also been 
investigated for these hybrid systems. [13, Recently 
zigzag graphene nanoribbons have been experimentally 
synthesized by selectively removing hydrogens from epi- 
taxial graphane. [soj Thus patterning offers a definite 
possibility of generating stable nanoribbons in the form 
of patterned graphene nanoribbons. 

Graphene being a potential candidate for spintronic 
device has led to a number of theoretical and exper- 
imental studies on the determination of the magnetic 
properties of graphene. [3l| The existence of a spin- 
polarized edge state has been predicted by a number of 
theoretical works, when graphene sheet is cut to have 
parallel zigzag edges (called zigzag graphene nanorib- 

bon or ZGNR). sElHii However, armchair graphene 

nanoribbon (AGNR) does not show such a spin-polarized 
state. The spin-polarized edge state in ZGNR is pre- 
dicted to be highly dependent on the edge geometry and 
may not be robust. [39[. Recently it was shown that 
transition metal (TM) termination at the edge trans- 
forms semiconducting ZGNR to a metallic one. It 
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was also shown that Fe terminated ZGNR shows an anti- 
ferromagnetic (AFM) couphng between the two edges of 
the nanoribbon. [4l[ Furthermore, it has been showed 
recently that the magnetic properties of transition metal 
doped ZGNRs are more robust than those moments aris- 
ing due to the edge geometry. j42| The most attractive 
feature of patterning graphene with partial hydrogena- 
tion is the similarity of creating bare carbon channels 
having similar edge state properties as those of graphene 
nanoribbons. Recently edge states of graphene/graphane 
interface has been studied using tight-binding approxi- 
mation and it was shown that edge state of an interface 
oriented along a zigzag direction enhances the effects re- 
lated to spin-orbit interaction, [i^ A recent study of 
zigzag graphene nanoribbon patterned on graphane using 
spin polarized ab initio calculations showed that the elec- 
tronic and magnetic properties of the graphene/graphane 
superlattice strongly depends on the degree of hydrogena- 
tion at the interfaces between the two materials. [44| 

Quite clearly the nature of electronic states in these 
two channels are different. The electrons in armchair 
channel are localized and non magnetic. However the 
electrons in zigzag channel are delocalized. A weak anti- 
ferromagnetic coupling across the edges exists which de- 
creases as the width increases. It will be interesting to 
investigate the effect of magnetic impurities in both these 
channels. Therefore we have carried out spin-polarized 
ab initio density functional theory calculations (DFT) to 
study the electronic and magnetic response properties of 
one and two magnetic impurities placed in these chan- 
nels. We have considered three different widths of the 
channels along with two to three possible Fe-Fe distances. 
We have also investigated the nature of the magnetic ex- 
change coupling between two Fe impurities. 

There is another aspect that needs attention namely 
the effect of spin orbit coupling. One may expect an in- 
crease in the orbital moment of nanostructures compared 
to that of the bulk due to reduced symmetry. An ex- 
perimental determination of magnetic moments includes 
both spin and orbital moments. Another point of rele- 
vance is the easy axis of magnetization. Therefore, we 
have calculated the orbital magnetic moments including 
spin-orbit coupling (SOC) in the Hamiltonian and have 
also discussed the strength of magnetic anisotropy along 
with the easy axis of magnetization. It is noteworthy 
that one of us has shown recently [45| that the magnetic 
anisotropy energy (MAE) of an Fe atom of Fe-porphyrin 
molecule deposited on a defected graphene lattice can be 
manipulated by straining the graphene lattice. In that 
case, the predicted MAE is between 0.4 and 2 meV. So, 
in the present case, it is interesting to study the effect of 
SOC induced changes in magnetism of the Fe impurities. 

The plan of the paper is as follows. In the next section 
we will describe in some detail our calculations followed 
by the results in Sec. IIIIl We will briefly recapitulate 
some of the important features of graphene/graphane in- 
terface. In Sub Sec. IIII Al we will discuss the effect of 
adsorbing a single Fe atom on both types of channels for 



different widths followed by the results on spin-orbit cou- 
pling induced orbital magnetic moments and magnetic 
anisotropy energies (MAEs) for a particular case in Sub 
Sec. IIII A II Then we discuss the results on the magnetic 
interactions between two Fe atoms in the channels. The 
conclusions are given in Sec. IIVI 

II. COMPUTATIONAL DETAILS 

We have used plane wave based density functional 
theory (DFT) as implemented in VASP code [M E^l 
for the calculations. The projector augmented wave 
(PAW) method [i^ within the generalized gradient 
approximation as given by Perdew, Burke and Ernzer- 
hof [13, IHH has been used for the exchange-correlation 
potential. The electronic wave functions were expanded 
using plane waves up to a kinetic energy of 450 eV. The 
spin-polarization approach was implemented in all cal- 
culations. For electron smearing, Fermi smearing was 
used and was kept fixed at 0.05 eV. The energy and the 
Hellman-Feynman force thresholds were kept at 10~^ eV 
and 0.005 eV/A respectively. We considered the chair 
conformer configuration of monolayer graphane, as de- 
scribed by Sofo et al . Q The pristine system was mod- 
eled using a supercell approach. The separation between 
the graphane layers in two consecutive unit cells in the 
perpendicular direction is kept as large as 15 A to avoid 
the interactions between the two graphane planes. For 
geometry optimization, a 5 x 5 x 1 fc grid was used 
whereas all x 11 x Ifc grid was used to calculate the 
Density of States (DOS) in the Monkhorst-Pack scheme. 
All atomic positions were allowed to relax where only 
the forces between the atoms are minimized and the vol- 
ume of the elemental cell is kept constant. This volume 
is equal to the volume of an equivalent graphane sys- 
tem. This is similar to the experimental situation where 
graphene nanoribbons are designed in a graphane sheet 
with fixed boundary. We have tested the influence of the 
cell relaxation and have found that the lattice parameter 
is decreased by less than 1%. which is insignificant for our 
present discussion. The use of a fixed boundary with the 
lattice parameter of graphane was justified in an earlier 
study [4J| , where it was argued that the relaxations of the 
internal coordinates is sufficient to describe the electronic 
structure and magnetism at the graphene/graphane in- 
terface. For the calculations of orbital moments and mag- 
netic anisotropics, the energy convergence was set to 10~^ 
eV. 

The binding energy Eh of Fe, is defined as 
Eb = [N * E{Fe) + E{Graphane)]~ E{FeN + Graphane) 
where 

- E{FeN + Graphane) is the electronic energy for 
equilibrated, optimized FeAr+graphane system. 

- E{Graphane) is the electronic energy for the opti- 
mized graphane geometry. 
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- E{Fe) is the electronic energy for the isolated Fe 
atom placed in a cubic supercell of length 10 A. 
Only the F point of the Brillouin zone was sampled 
in this case. 

- A'^ is the number of Fe atoms. 

III. RESULTS AND DISCUSSIONS 




(a) "armchair" channel 




(b) "zigzag" channel 



FIG. 1. (Color online) Decoration of hydrogen (a) along the 
diagonal of the unit cell leading to an "armchair" channel 
and (b) along the edge of the unit cell leading to a "zigzag" 
channel. In the figure, yellow (light shaded in print) balls 
are bare carbon atoms, turquoise (dark shades in print) balls 
are hydrogenated carbon atoms and red (small dark in print) 
balls are H atoms. For both armchair and zigzag channels, 
the nanoribbons are 3-rows wide. 

We briefly recapitulate the interesting result of two 
types of graphene/graphane channels. These two chan- 
nels are schematically shown in Fig. [T] Figures [H^a) 
and[ljb) show "armchair" and "zigzag" channels created 
by removing hydrogen atoms along the diagonal and the 
edges of the supercell respectively for a specific width of 
3 rows. While the "armchair channel" is non magnetic 



TABLE I. Width dependent energies and magnetic moments 
of the channel systems with a single Fe impurity. Ei,, /itotai 
and /ipc denote the binding energy of Fe, total magnetic 
moment of the system and magnetic moment on Fe site 
respectively. 



Single Fe atom 


Channel Type 


Channel Width 


Eb 
(eV) 


/itotal 

(/is) 


/iPo 

(/is) 


Armchair 


3-rows 


1.6 


2.0 


1.99 


5-rows 


1.10 


2.06 


2.10 


7-rows 


0.80 


2.00 


1.98 


Zigzag 


3-rows 


1.38 


2.1 


2.5 


5-rows 


1.42 


2.01 


2.55 


7-rows 


1.45 


2.01 


2.56 



and has a width dependent bandgap, the "zigzag chan- 
nel" interfaces are magnetic. The lowest energy structure 
is antiferromagnetically coupled across the edge but as 
the width increases, the FM and AFM coupled edges be- 
come almost degenerate e.g., FM-AFM energy difference 
decreases to 3.63 mcV/edge C atom from 6.57 mcV/edge 
C atom as the width of the zigzag channel increases from 
3-rows to 7-rows. Beyond a certain width, the exchange 
coupling between the edges becomes zero with C mo- 
ments along each edge ferromagnetically aligned. 

A. Single Fe atom on bare channels 

We begin by presenting the results of a single Fe atom 
placed on three different widths of these channels. It is 
known that for pristine graphene, the energetically favor- 
able position of an Fe atom is at the hollow site of the 
hexagonal ring with an Fe to carbon distance of 2.1 A. 
[52h54| Our calculations show that the minimum energy 
position of Fe is dependent on the type of channel. For 
the armchair channel the lowest energy position of Fe is 
at the hollow site equidistant from the interface for all 
the widths. However in the zigzag channel, the mini- 
mum energy position of Fe is when it is at the hollow site 
of bare carbon atoms nearest to the interface. This is 
because of the anti-ferromagnetic nature of the edge cou- 
pling. For the 3-rows channel, the central row is nearest 
to both the interfaces and Fe atom prefers to sit there. 
However as the width increases, the Fe atom shifts to- 
wards the interface. The binding energy of the structure 
where Fe is in the interface is higher by 0.31 eV/Fe atom 
and 0.36 eV/Fe atom respectively for the 5-rows and the 
7-rows channel than the Fe atom at the center of the 
channel. The position of Fe is slightly asymmetric w.r.t. 
the surrounding hexagon due to the stretching of under- 
lying carbon bonds and the Fe to C distance varies from 
2.06 A to 2.18 A. 

Now we will discuss the energetics and electronic struc- 
ture of a single Fe atom placed in both the channels for 
all widths. Table U shows the energies and magnetic mo- 




FIG. 2. (Color online) Spin polarized DOS for single Fe atom 
placed on the 3-rows "armchair" channel, (a) Total DOS with 
(blue solid line) and without (red dashed line) Fe. (b) Site 
and angular momentum projected DOS for Fe surrounding 
six C atoms. The orange (in print, light grey) shaded region 
shows the pz component, (c) Site and angular momentum 
projected DOS {d components only). 



mcnts of the Fe atom placed in these channel systems. 
From the table, it can be seen that as the channel width 
increases, the binding energy of Fe atom for the armchair 
channel decreases while for the zigzag channel it remains 
more or less constant. We also note that the total mag- 
netic moments for both the channels for all the widths 
are 2.0 fj,B, consistent with the value for magnetic mo- 
ment of a single Fc atom on graphene. [53| However, 
the on-site local magnetic moment on the Fe atom is 
0.5 fiB higher in the zigzag channel (see the discussion 
next). It is also observed that the binding energy of Fc 
in all three different width zigzag channel structures are 
about 0.2 eV higher than the same in pure graphene. [s^l 
Thus the mixed sp^ — sp^ character increases the binding 
energy of Fe in the zigzag channels. 

We further analyze the result by examining the DOS 
and spin densities of both the channels. In Fig. [2] we 
show spin polarized total DOS with (blue solid line) and 



FIG. 4. (Color online) Spin density plot of a single Fe atom 
on the 3-rows "armchair" channel. Red (in print dark shade) 
is positive spin density and yellow (in print light shade) is 
negative spin density. 

without (red dashed line) Fe (fig. [2ja)), site and angular 
momentum projected DOS for Fe and Fe near-neighbor 
C atoms for 3-rows armchair channel (fig.^b-c)). As can 
be seen from the DOS, there are induced states due to 
spin-down electrons appearing in the gap below the Fermi 
energy. These states are due to Fe d spin-down electrons. 
As a consequence, the gap reduces significantly. This gap 
is in between spin- up and spin-down electrons. Analy- 
sis of total DOS for 5-rows and 7-rows armchair channel 
(shown in fig. [3]) indicates that this feature is also present 
as the width of the channel increases from 3 to 7 rows. 
However the value of the gap depends on the width of 
the channel. 

The nature of spin densities (shown in fig. |4]) can be 
analyzed with the help of site and angular momentum de- 
composed density of states (henceforth called as PDOS) 
(shown in fig. [DJb-c)). The spin density plot shows that 
the majority of spin-up density is localized only on Fe 
atom and Fe is seen to induce the spin-down density on 
the surrounding six C atoms. There is a large localized 
magnetic moment of the order of 2 /x^ on Fc (see table H]) 
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and a very small negative magnetic moment on the sur- 
rounding six C atoms -0.02 fiB/atom). Interestingly 
there is an induced up spin density in between the car- 
bon atoms of Fe surrounding ring. Analysis of PDOS 
shows that the Fe d^z and dyz of up spin channel inter- 
act with the TT bonded system formed by pz orbitals of 
the surrounding carbon atoms (shown in fig.[31Jb)). This 
induces polarization on the ring so as to generate neg- 
ative moments on the C sites and positive moment in 
between the C atoms (see spin density in fig. 2]). It is 
to be noted that the magnetic moment of Fe is reduced 
from for isolated atom to 2fiB for the system. The 
analysis of angular momentum decomposed DOS of Fe 
(fig. [21(c)) shows a strong hybridization of Fe 3d states 
with the Pz orbitals of the C atoms. The DOS shows that 
all five spin-up 3c? states and three spin-down 3d states of 
Fe are occupied. Both spin-up and spin-down 4s states 
are above the Fermi energy indicating that two electrons 
are transferred from atomic Fe 4s states to spin-down 
3d states of Fe-channel system. This causes the reduc- 
tion of magnetic moment on Fe site which is consistent 
with the earlier result, [s^. For larger widths ( 5-rows 
and 7-rows), we have also analyzed the site and angular 
momentum projected DOS (figure not shown) and the 
behavior of the spin densities are found to be similar in 
nature. 

Fig.[S]-(I) shows spin polarized total and site projected 
DOS for the 3-rows zigzag channel. The red dashed line 
in the total DOS is for a bare zigzag channel. In contrast 
to the armchair channel, the DOS for the edge channel 
(fig- EK^-)) shows a finite DOS near the Fermi energy. 
This is mainly because of spin-down electrons from the 
Fe site and bare C sites with peaks below and above 
the Fermi energy. The detailed examination of site and 
angular momentum projected DOS shows that the down 
spin peak just below the Fermi energy comes from both 
Fe dz2-r^ orbitals and the Pz orbitals of the interface bare 
carbon atoms. The contribution in up spin peak just 
below the Fermi energy is from delocalized Pz orbitals of 
edge bare C atoms as well as a small contribution coming 
from TT bonded C atoms surrounding the Fe atom. 

Our analysis also show that the effect of magnetic 
impurity is felt up to 4th nearest neighbor interface C 
atoms. It reduces the on-site magnetic moments on 
edge C atoms and the maximum reduction is ~ 15% on 
the nearest site. This is a sharp contrast compared to 
the armchair channel where the Fe impurity affects only 
up to 1st nearest neighbor. Another interesting differ- 
ence is in the induced spin density on the hexagonal C 
ring nearest to Fe. In the present case this ring is not 
anti-ferromagnetically ordered. Only three C atoms (be- 
longing to same sub lattice) show significant down spin 
density (see in fig. [5]-(II)). This is understandable from 
the fact that the delocalized electron density present in 
the channel interacts with the Fe atom. It reduces the 
amount of charge transfer from Fe 4s states to Fe 3d 
states, resulting in an increment of local magnetic mo- 
ment on Fe to 2.5/iB. This also leads to a somewhat long 




(11) 

FIG. 5. (Color online) (I) Spin polarized DOS for a single 
Fe atom placed on the 3-rows "zigzag" channel, (a) Total 
DOS with (blue solid line) and without (red dashed line) Fe. 
(b) Site and angular momentum projected DOS for Fe sur- 
rounding six C atoms. The orange (in print light grey) shaded 
region shows the Pz component, (c) Site and angular momen- 
tum projected DOS (d components only). (II) Spin density 
plot of a single Fe atom on the 3-rows "zigzag" channel. Red 
(in print dark shade) is positive spin density and yellow (in 
print light shade) is negative spin density. 
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ranged perturbation in the underlying C lattice. 



1. Role of spin- orbit coupling 

As discussed in the introduction, we have included 
SOC in the Hamiltonian to calculate the orbital moments 
and magnetic anisotropy energies for the case of 3-rows 
wide armchair and zigzag channels with a single Fe atom. 
The calculated orbital moment of the Fe atom in the arm- 
chair channel is 0.09 ^s, which is enhanced compared to 
the calculated value for bulk Fe 0.05 jU without or- 
bital polarization) in the bcc phase. This is expected due 
to the lowering of symmetry in the present case. How- 
ever, the enhancement is not dramatic due to the strong 
crystal field produced by the C atoms in the hexagon 
to which Fe is chemically bonded. The calculated easy 
axis of magnetization is in-plane with MAE equal to 19 
meV/Fe. This strong magnetic anisotropy should hold 
the Fe moment in the plane of the substrate. The or- 
bital moment along the 001 direction (hard axis) is 0.02 
fiBi consistent with the fact that the orbital moment is 
reduced along the hard axis. We have observed an in- 
teresting effect in the case of the zigzag channel. In the 
presence of SOC, some of the C atoms have spin moments 
of the order of 0.1 /is. Therefore the total moment of the 
unit cell in much increased (4.5 ^lb) compared to the case 
without SOC. Although the Fe spin moment stays more 
or less the same, the additive contribution from C atoms 
gives rise to this enhancement of the total moment. Also, 
in the case of zigzag channel, the MAE is smaller 4 
meV/Fc) with an in-plane easy axis. The in-plane or- 
bital moment amounts to 0.1 ^b, which is slightly larger 
than the out-of-plane moment (0.08 /is) The topic on 
the magnetic anisotropy of adatoms on graphene and the 
tuning of MAE with an electric field has been discussed 
very recently. [55| 

B. Magnetic Interaction between two Fe atoms 

Now we will discuss the case of two interacting Fe im- 
purities. As expected, it turns out that the response of 
the channel electrons are different for the armchair and 
zigzag channels. Table ITTl shows the binding energy of Fe 
atoms, total magnetic moment, on-site magnetic moment 
on Fe atoms and exchange energies (wherever possible). 
For armchair channel we have considered two possible Fe- 
Fc distances for all the channel widths. In this channel 
the two magnetic impurities (Fe) are weakly interacting 
as evident from the weak variation in the binding en- 
ergy of Fe atoms as a function of distance coupled with 
the fact that the total magnetic moment remains almost 
unchanged. Consequently the exchange energy {E*^^) is 
also negligible. The spin densities are also similar to the 
isolated single Fe atom in the armchair channel (additive 
in nature). Our calculations show that these features are 
also present for the larger widths. 



However, the examination of DOS shows some inter- 
esting features. The DOS of two Fe atoms in the 3-rows 
armchair channel (fig.El^a) and fig.ElJb)) shows that both 
Fe atoms induce states below and above the Fermi energy. 
It is interesting to note that the DOS shows a gap just 
below the Fermi energy in the spin-up channel and just 
above the Fermi energy in the spin-down channel. The 
observation of DOS for small Fe-Fe distance (see fig.lH^a)) 
suggests the possibility of a spin gapless semiconductor 
material (SGS). It can be seen that there are energy gaps 
of 0.56 eV and 0.45 eV in the up and down spin channels 
respectively. However, the cross spin gap between down 
spin channel and up spin channel is very small. There- 
fore, at the Fermi energy, an electron has to flip its spin 
to go to the conduction band from the valence band. The 
resulting transmission is 100% spin polarized. Recently 
Wang has discussed such a possibility in SGSs. [5^ It 
was further noted that SGS can be obtained by doping 
magnetic impurities to semiconductor. In the present 
case, a similar situation occurs. As the distance between 
two Fe atom increases, this feature can also be seen there 
(see fig. m^b)). However, as the channel width increases, 
this feature weakens. Thus there is a possibility of SGS 
for a narrow strip of armchair channel. 

In contrast to the armchair channel, the binding energy 
as well as the magnetic moment are distance dependent 
in the zigzag channel. For the smallest Fe-Fe distance 
for all the widths, the two Fe atoms forms a dimer. It is 
evident from high values of the binding energies as well as 
the total magnetic moment of the system (see table |ll]) . 



TABLE II. Width dependent energies and magnetic moments 
of the channel systems with a pair of Fe atoms with possible 
Fe-Fe distances. The distances between two Fe atoms (d) are 
in Angstrom, binding energies of Fe in electron volt {Eb). 
Total magnetic moment of the systems (/^totai) and magnetic 
moment on the Fe atoms (/jfc) are in the units of Bohr- 
magneton and the exchange energies {E""^ = — E^^'^') 
are in electron volt 



Two Fe atoms in the channel 


Type 


Width 


(A) 


Eb 
(eV) 


Mtotal 


(/is) 


(eV) 


armchair 


3-rows 


4.38 


3.13 


4.02 


1.98/2.01 


0.009 


8.77 


3.18 


4.01 


1.99/1.99 


0.001 


5- rows 


4.38 


2.25 


4.11 


2.10/2.10 


0.018 


8.77 


2.19 


4.08 


2.11/2.12 


0.008 


7-rows 


4.38 


1.75 


4.06 


2.04/2.03 


0.014 


8.77 


1.64 


4.01 


2.01/2.01 


0.016 


zigzag 


3-rows 


2.17 


4.12 


5.96 


2.98/2.98 


0.457 


5.08 


2.67 


7.0 


2.59/2.6 


0.058 


7.60 


2.67 


4.26 


2.42/2.43 


0.038 


5-rows 


2.17 


4.09 


5.89 


3.01/3.01 


0.834 


5.08 


2.73 


4.14 


2.53/2.53 


NA 


7.60 


2.80 


4.06 


2.51/2.51 


NA 


7-rows 


2.17 


4.10 


5.82 


3.01/3.01 


0.795 


5.08 


2.80 


4.07 


2.53/2.54 


NA 


7.60 


2.87 


4.02 


2.52/2.52 


NA 



7 



For comparison we note that the total magnetic moment 
of Fe-Fe dimer is 5.94 fj,B- The on-site magnetic moments 
on Fc atoms are also very high ('^ 3 /is on each Fe atom). 

Table |lT] also indicates that there is a significant varia- 
tion in total magnetic moment for the Fc-Fc distance of 
5.08 A for all the widths. Although the individual mag- 
netic moment on Fe atoms are similar and of the order 
of 2.5 Hb, the total magnetic moment reduces from 7.0 
fiB (3-rows) to 4.14 and 4.07 yUs for the 5-rows and 
the 7-rows channels respectively. In this case the spe- 
cific position of 2nd Fe atom induces negative spin den- 
sities on the surrounding C atoms which in turns affect 
the spin densities of both edge row C atoms (see fig.[71). 
As a consequence, the channel becomes ferromagnetically 
coupled across the edges and develops a large total mag- 
netic moment. For the larger widths, as the Fe atoms are 
placed nearest to one edge, the interaction between the Fe 
atoms and the other edge is much weakened. Therefore 
the edges remain weakly anti-ferromagnetically coupled 




-3-2-10123 
E - E, (eV) 

(a) Fe2 DOS - 3-rows (4.38A) 




-3-2-10123 



E - E, (eV) 
(b) Fe2 DOS - 3-rows (8.77A) 

FIG. 6. (Color online) Spin polarized total DOS plots for 2 
Fe atoms in the 3-rows armchair channel. We've plotted the 
graphs within -3 eV to 3 eV to bring out the feature of DOS 
near the Fermi level. The numbers in the parenthesis are the 
distances between the two Fe atoms. 




(a) Fe2 Spin 3-rows (2.17A) 




(b) Fe2 Spin 3-rows (5.08A) 




(c) Fe2 Spin 3-rows (7.60A) 



FIG. 7. (Color online) Spin density plots. Red (in print dark 
shade) is positive spin density and yellow (in print light shade) 
is negative spin density. 



reducing the total magnetic moment. 

As expected, the exchange energy E'^^ is stronger for 
the smallest Fe-Fe distance (2.17A) for all the widths 
(see table. |ll|. For 3-rows width, as the Fe-Fe distance 
increases, the exchange energy value decreases to 0.038 
eV. Thus there is a stable ferromagnetic state for all the 
Fe-Fe distance with possible transition temperature ^ 
350 K. However it was not possible to estimate E"^^ in 
other widths because we could not get properly converged 
anti-ferromagnetically state for longer Fe-Fe distances. 

From the above discussions, it is evident that for the 
zigzag channel, the magnetic atoms interact strongly. 
Thus one can envisage to achieve a magnetic graphane 
lattice by depositing magnetic atoms (Fe in the present 
case) in the metallic channel of bare carbon atoms. In- 
deed our preliminary calculations shows that it is possi- 
ble to obtain stable linear Fe chain by placing Fe atoms 
nearest to the interface of the zigzag channel. The total 
magnetic moment and the binding energy was found to 
be 2.42 /is/Fe atom and 2.14 eV/Fe atom respectively. 
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It may be possible to form novel magnetic nanostructures 
(chains, islands etc.) with the help of scanning tunneling 
microscopic tips and study magnetism as a function of 
the width of the bare carbon channels. Therefore, an ul- 
trathin device with novel functionalities may be achieved. 

IV. CONCLUSION 

In conclusion, our detailed density functional investi- 
gations have revealed some novel electronic and magnetic 
response of single and a pair of Fe atoms placed in pat- 
terned dehydrogenated channels in graphane. Our work 
shows that the response of the two channels - the "arm- 
chair" and the "zigzag" channels are different. In the 
armchair channel, the magnetic response of the single Fe 
atom is localized in nature while for the zigzag channel it 
is fairly long ranged. The easy axis of magnetization is in 
the plane of graphene with a large magnetic anisotropy 
energy of 19 meV/Fe atom in the case of the armchair 
channel. As compared to armchair channel, the binding 
energies and the total magnetic moments with two Fe 
impurities are distance dependent in the zigzag channel 
suggesting that the pair of Fe atoms are rather strongly 
interacting in this channel. The magnetic coupling be- 
tween a pair of Fe atoms is very weakly ferromagnetic 
in the semiconducting armchair channel whereas a rela- 
tively stronger ferromagnetic coupling is observed for the 



closest separation between the Fe atoms in the delocal- 
ized zigzag channel. Our result for the armchair chan- 
nel with two Fe atoms shows a possibility of having a 
novel spin gapless semiconductor material by doping ap- 
propriate magnetic impurities in partially hydrogenated 
graphene of a particular width. The magnetic interac- 
tion between the impurities in partially hydrogenated 
graphene with conducting channels shows an oscillating 
ferromagnetic and anti-ferromagnetic coupling across the 
edge of the channel for a particular width of zigzag chan- 
nel. The possibility of realizing an ultrathin device hav- 
ing novel properties is discussed. 
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